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A  STUDY  CF   FPiiSCUllL   GCUTROL   FOR    ^  BLOWDOWN  VIND    T^JNl^Ii,L 

The    proble/a   involved    in    this  investic^^tion  was 
tc  detertiiinb   tiie  characteristics  of   a  syateia  that  v.oula 
effectively   control   the    settling  chaiubor  pressure  within 
Cfcrttiin  arbitrary  limits. 

The   study   wais    coiicferntid   only  vdth   one   particu- 
lar   installation.     Hov/evur,    the  methods  eni ployed   should 
be  applicable    to    ony   si.iilar   problen. 

As    u  rcisalt    of    tnis   investiii^a  tion   it  ..:.ay    by 
concluded   thi^t   x^i'^ssure   control   to  v;ithin  one   pound  per 
square    inch  :;iay    be    obtained  v>^ithin  apDroxirately   one 
second.      This   assa;:\t}o    one   Uiit   oj.    '■    otfrvo^i',:;cxiciuii..i  in 
which  the   control   valve   is  moved   by  means  of  a  servo- 
motor  fictinr   t''Touf>h  a   gear    train  and   rac.^:.      The 
[/^Li^rul    realoi'yfc;    Ox     the    serV0;'it;Ciiaiiic.:LL  f.i'^   xiirttd    in 
Table  II.     A  notor   of  such  characteristics  as  are 
required  to  provide  this  control   is  described   in  this 
i'  ep  0  rt « 

It   is  recoiin.'iended   that   such  a   servomechaniSia 
be   desi{^ned   end    installed    in   the   subject  wind  tunnel 


in   the   place    of    the   present   installatioi..      It    is    ex- 
pected  that   this  will    be   accojiplished   in  the  near 
future. 

This    investication  was    conducted    in   the 
laboratories  of   the   Aeronautical  iin^ineerin^;  Depart- 
iTient  of    the  University    of  iTichi^an  by  F.N.   Kanaga 
and  D.B.   Rice.      It   is    sub. ratted  as  a   thesis    for 
credit    towards  the  degree  of  Master  of   Science. - 
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HITRCDUGTION 

Tile  problei.i  to  be  investigated  is   the  control 
of   the    air   pressure    in   the   settling  chamber    of  a   blov.- 
dovn  wind   tunnel.      Thy   particiiiar  wind    tunnt;!   to   be 
studied   is  already   iiistalled   in  the  East  Engineering 
Building:  of   tlie   University   of  Michii:ian  at   Ami  Arbor, 
Michic^n.      A  photo£;rfcip]:i  of    the   installation   is   shown 
in  Figure  1,   and   soriie   of    the   physical  measurements 
are    in  Table    I, 

The   principal   components  of   the   wind  tuniiel 
are   the    storage   tank,    the   driers,    the  quick-opening 
valve,    the   control   valve,    tht;   settling  chamber,    the 
nozzle    and    diffuser,    una  the   control   system.      Air 
pressure   at  about  95   pounds   per    s^itiare   inch  ga£-e    is 
furnished   by   the  building:   air    supply   system. 

The    perfect  control    systeja  in   tiiis   case  would 

» 
be    one    thb  t  v/ould  ueasui'e   the    instantaneous   air  pres- 
sure   in   the   settling:;   chamber,   coapare   it  with    the   de- 
sii'er    pr'essure,    and  then   actuate    the  control   valve    so 
that  the  pressure  wo  aid    be   changed   to   exactly   tnu 
desired  pressure   in  an   arbitrarily   short  period  of 
tijie.     When   the  v/ind  tunnel   is  employed  for   sun-r-^>onic 
laboratory  tests,    the  length   of  blov/dov/n  runs  v^ry 
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fr'J-i  about   ;>,;■    to   iuuout   25    seconds.      Tri3  vuriation 
depends  primarily   upon   the  value  of   the   desirfccl   pres- 
sure   in  the    settling   chainbfer,    providing    the    initial 
supply   tank   pressure    is   constant.     To  make   inaxiiauiii 
usf   of    this    short  tiine,    it   is  desired  that   the  pres- 
sure be    stabilized  within  one   second.     Frorn  exper- 
ience,   it  v;a3   decided   that   control   of   the   settling; 
tank  pressure  to  v/ithin  plus  or 'minus  one  pound  per 
square    inch  would   be   as   close  as  vwnuld   conceivably 
be  required    in  tiie   laboratory.      Therefore,    the. 
arbitrary  goal    set   for  the  control   systeia   is  to  re- 
duce   the  pressure    error    (instantaneous  settling 
cha:iber   pressure  lainus   desired  pressure)    to  within 
plus  or  Liinus  one  pound  per    square    inch  within   one 
second • 

The   present  control    systej-i  installation  is 
of   a  type   in  general  use   in  the   chemical   industry  for 
tej-uperature  control   and  was  modified  for   this  pressure 
control   application.     This   installation   is  not   satis- 
factory  for    this  wind  tunnel  lor    several  reasons.     The 
primary    reason  is   that  the   pneumatic  portion  of  the 
circuit   is   such   that   the   cohtrol    air  pressures   are 
too   small   to   move  the  diaphragm  actuated  control- valve 
vvith  sufficient  velocity  to   ^ive  close  control  of   the 
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sett  lino;   chamber   pressure,      witri   the   control  system 
on  "automatic",    it  re^iUired   20.7   seconds   for   the 
valve   to    open   fully,    and  13.5    seconds  for   it   to 
close.      Increasiiig  the   control   air  pressures  would 
involve  very   extensive  rev^orkini?;   of   the  pneuaiiatic 
system  and    does   not  appear    to  offer    a   substantial 
imj^rovement  due   to   the    prtdoui-y  1  Lai  tat  ions  of  the 
valve  diaphragm. 

In  view  of'  the   above,    another  metnod  of 
pressure    control  hc.s  been  investigated   and  made  the 
subject   of   this  report.      The    e^mployment   of   an  electric 
iiotor    to   drive    the  valve    stem   througti    a  e^ear    train   and 
rack  was  dec-ided  upon   as  a  likely   possibility    for   pro- 
viding  valve  movement,    this  machinery    to   be   the  load 
for    some  sort   of   servo uiechan ism.      ^  sketch   of  this  pro- 
posed  arrangement    is  shown  in  Figure  2o 

The  method  of   investigation  v.as  to  determine 
the   forces   that   act  on   the  valve    end    of   the  valve   stem, 
and   then  to  determine  the  control  forces  necessary   to 
overcome  these  valve  forces  and   to   move  the  valve   to 
a  position   corresponding  to    tiie  desired  pressure   in 
a   stable   manner  vdUxin   the   desired   tiiAe  limit. 

The  major    tool   employed  in  this  study   vvus 
the    electronic   differential   analyser,    (electronic 


cnaloG   computer)   vvjl-':    win    >'^   referred   to  hereafter 
in  this  report   as    the   "coiiouter "..      Inforiaation  con- 
cerning  this   coaiputer  iiiay    be  found   in  references 
(-)  .    (  .^)  ,    (c)  ,    (d)  ,     1^4    (e)  . 

.■%  differential  ©'-•uation  was  developed  which 
eaiLouits    tlie   forces  that    act   on   tiie  valve.      By  iaaking 
use   of  the    various  functions  of   the   operational   ampli- 
fiers,   which  constitute   the   cauputer,   for    sign   ch ending, 
nultiplication  by  a  const^-int  factor,    integration  and 
summation,    tlie  differential   equation  was   duplicated  on 
the    computer.      The  voltages   at   various  points  witiiin 
the   coiaputer  represent    tJie  values   of    corresponding 
functions  within   the    equation,    and    these   voltages   are 
fed  through),   a  Brush  direct   current  ajaplifier   to   a 
Brusii  jTia£;netic  recording   oscillograpii,    or   recorder, 
so   that  the   variation  with  tii-ie   of  the  various  funct- 
ions   in  the  ecuation  wbs  recorded  on   tape. 

This   investigation   is  limited  to   one   specific 
wind  tunnel.      The  methods  used,    however,    siiould  be 
applicable   to    other    such   systems.  *■ 

This  study  was  .aade   during   the   1949-1950 
school  year   in   the   laboratories   of    the   Aeronautical 
Engineering   Depai'tment   of   the  University   of  Michi£^an, 
The    autliors  wish    to  acicnowled^e    and   express  appreciation 
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for   the    invaluable  guidajice  and   sug^^^stions   eiiven  by 
Dr.   L.L.    Rauch    of   the   Aeronautical   Ln^ineering  Depart- 
iTie  nt . 

PHOCEDURK 

wniD   TUTTII'CL  DATA 

The   v^ii'lation  of   settling  chamber   pressure 
with   valve   posit iori   as  shown   in  Plate   I  was  found  by 
presetting  the  control  valve   by  laanual  laeans    to  various 
positioiiS   and   then  reading    the  inaximujii  settling  ciianiber 
pressure  reached  \tien   the    quick-openin^j  valve  was  opened. 
These  values  were   read  from  a   servo-type  electrical   in- 
dicator Y/hich  was  actuated   by   a  differential    transfonaer 
attttohed  to  a  tourdon  tube  located  on  the  settling  cham- 
ber.     The    initial   tank  pressui^e  was  constant   for   each 
set   of  runs. 

The    coulorab  or  dry   friction  and  the  weight  of 
the  movinc  parts  of   the  valve  were  determined   as   indicat- 
ed  in  Figure  3.      The  forces  Dieasured  with   the   valve   barely 
jaoving  up  vjere  weii^ht   plus   friction,    and  with  downward 
notion    tiiey  were  v;eigiit  minus   friction.     From  this  the 
weight  was  determined   to    be  10, 3S  pounds   and   the   cou- 
lomb friction  to  be  1,23  pounds. 
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Lift  forces   on  tne   valve    due    to   air  passing 
through   the   valve  were  measured  by  usin£^   apparatus  set 
up   as   in  Figure  4,      a  plot   of  lift   forces  versus  valve 
position  is  on  Plate   II. 

No  Qtteapt  was  aade   to   measure   the  effect   of 
viscous  friction  of  the  valve.      The  force   necessary 
to    start   the  valve  moving;   very  slowly  was  about   one 
pound  ci't^ater   than   tlie   coulomb   friction  force.      This 
force  difference   is  sometLoies    called   "sticiction" . 

CCr/IPhNS/iTING   FOR    ffiE  LIFT  FGHCE 

Referring  to  Plate   II,    it  will   be  noted  that 
the  lift    force   due    to    tXx&   action  of   the   air  passin^^ 
throut^h  the   control   valve  has   been   approximated   as   a 
linearized   function  of  valve  position.      The  weight 
of  the  valvfc!   and   stera  will   cancel   out  10, 3S  pounds 
of    this   force.     However,    this    still  leaves  a  lar^e 
force    to   be    supplied   by   tixe    control    system  when  the 
valve    is   being  operated  near    the   fully  opened  position. 
By   addint;   a   six  pound  weight    to  tiie  valve    stem,    the 
lift  force    is    btUanced   by   v/eicht  v^hen   the  valve   is 
one    inch,   or   two-thirds  of   the  v;ay,   open.     This  six 
pound  weight   is  shown  in  Figure   2  us   an  addition  to 
the   weight  of    tiie  valve    stem.      It  v^ili   be   noted   that 
\^h.en  the   control  valve   is  fully  opened  the   uncompen- 
sated lift   force  will  be   six  pounds. 
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SELECTION  OF  MOTOR 


From  Computation  I,   vjhich  is  based   in  part   on 
a  rough  assiimption  of  accelerations,    horsepovrer  reouire- 
iaent3  for  a  rnotor  to  dr^. -•■   •*  atrol  valve  v;  ere  found 

to  be  0.000564  H?.      A   servomotor  thiit  vrould   deliver  about 
10  times  this  horsepov/er  at  riiaximijjii  output  was  tentative- 
ly  selected.      This  --      tho   "jiehl,   I'svy  type  C3A   211052 
motor  for  v/hich   the   characteristic   curves  are  shovjn  on 
Plate   III.      If  geared  30  that  the  motor  would  accelerate 
to   2000  rpLi  in  1/2    second,    .i ..    ras  found  that  an  average 
force  of  20,5     pounds  i-/ould  be   exerted  on  the  valve  stem, 
varying  from  23.^  pounds  at   aero   rpm  to  15.3    pounds  at 
2000  rpLi.      It  was   estimated  that  the  forces   involved 
would  be  below  these  figures.      The  gear  ratio  necessary 
to  obtain  the  average  force   is  69.7  to  1,   that   is,    the 
gear   that  drives   the   rack   should  turn  one  revolution 
v/hile  the  motor  turns  69.7  revolutions. 

After  a   solution  to  the  nrobleri  was  obtained, 
it   vao  apparent   that  --     :'-^tor  wliich   CDuld   exert    a   greater 
torque  vrould  be  desii*eable  in  the   interests  of  a   larger 
nargin  of   safety.      The  Diehl  Navy  type  GDA  211426  or 
211454,   which  v/as  found  t  ■  ^^-  '.r^-   -bout   doubl'--   ^^■- ^   p--^'-'^v 
output  of  t^ie  original  motor,    v/o.s   chosen.      Character- 
istic curves   for  this  motor  are  shovm  on  Plate   IX. 


TI-Ii:   DIFFERENTIAL  E--JJATICN 

Tlie   forces  on   the   valve   stem  were   aumiaod  up 
in  a  differential    equation  as  follows: 


^'ifl  ^^^'^  ^f-  -^  ^^P-  ^  ^^£?e^^  +  K.Pd  =  0 


ni?^-#-L.  ^+c..  x  +  r\.D_+n.n    4.K    l»^<^-^J_u'n._n  ^)-, 


where: 

X  is  valve  position  in  feet  ueasured 

« 

from  closed    to   open.     Full   valve   travel 

is   1.5   inches 'for    tiiis  valve.     For   the 

coordinate   systeia,    see  Figure   5. 
^  is   the    iiiass   in  slugs  of    the  moveable 

valve  parts  plus  the    six  pound  v;ei{;ht 

previously  mentioned. 
2  is  the   instantaneous  pressure   in  pounds 

per   square  inch   in  the    settling   chamber, 
P  is    the    desired  pressure   in  pounds  per 

sciuare   inch   in  the   settlin/^   chamber. 
Pe  ^^^  the   pressure   error   in  pounds   per 

square  inch   in  the   settling  chamber, 
cr    {   p  -   p   ). 
£x  is   thr-   coulomb  friction  force   in  pounds. 

X  is  a   unity  term  whose   sign   furnishes  the 

direction  of    the  coulomb   friction  force. 
Cp  is  a  coefficient   of   viscous  friction, 

with  units  of  pound   seconds/feet. 


—1'  —2  ^^^  — ?     ^^^   conotiints  v/hich  when  multiplied 

by  Pe»  Pp  ^^cl    I  p  dt  respectively 
represent  the  forces  in  pounds  to 
be  supplied  to  the  control  valve 
stem  by  the  control  Liechanism. 
K,P,  is  a  disturbing  force  in  pounds  on 

the  valve  itself. 
In  order  to  reduce  equation  (1)  to  a  function  of 
one  variable,  the  relation  obtained  frora  Plate  I  is  used: 

(2) 
p  =  n.  ■*-  n3,x 

v/here  n-,  is  the  ordinate  at  x  =?  0  and  n^  is  the  slope  of 
the  line.  ■(   ^^.ef erring  to  Plate  I,  iu  is  seen  that  equation 
(2)  is  valid  only  for  the  portion  of  the  carve  belov;  a 
valve  opening  of  about  0,575  inches.  The  nonnal  operat- 
ing region  for  the  control  valve  is  in  the  strain;ht  line 
portion  of  the  curve,  therefore  the  linear  relationship 
between  settlinf^  chanber  pressure  and  valve  opening  was 
used,   -^hus  from  the  above: 

P..-P-P  (3) 


p,  =p  =  n,x  ,  i^^pe/Hx 

3ubstitution   of    (4)    and    (5)    into    (1)    give; 
differential  equation  v/ith  p^/n^,   as  the  variable. 


(5) 


By  (iefinint;  K.  n     as  k   ,   K  u     as  k     and  K  n     us  k   , 
12  1<,2  2  32  3 

eiiuytion    (6)    becoxiies  nore   siiniUe    in  fovm: 


^''  '^..  V ;  *■■    .  ^  ■  (7) 


The   v8rioas  terms  or    the   differential   equation  represent 
physical   effects   as   followa; 


'^(PqX  is    tiib   effect   of  the  imas  of   tne  laoving 

'-  parts  of   the    control   valve   i^lus   tiie    six 

pound  compensatinti  weight, 

C-    6  is  the   coulo/nb  friction  force,    of  magni- 

1  2^  '  *  ^ 

•Pgl  tude  C,    and   vc^^^  L.i.-   in  sign  according];  to 


the    sign  of   d   • 

"e       . 

C^/p   V  is  thfc   viscous   friction   force    caused   by 

'2'  the  motion  of   the  control   valve. 

is    the   derivative   term  or    the   error-rate 
viscous   damping  force   to   be   supplied   to 
the   valvt^   stexa  by   the    controller. 
k    /p  I  is  the   direct  proportional  force   to   be 

^^2^  supplied  by   the    controller. 


'a 


'-'\& 


lat  is    tfje    inte,jral    control   force   to   be 

Q  "2  supplied  by   the   controller. 

_   1  ?   _    f 


K. ?^        is  a  disturbing  force  applied  to  the  valve 

stem  itself,  and  may  re'^rcsent  the  unbalanced 
lift  force  on  i,i.e  vir^lvy  v/hicii  ".a^  previously 
discussed  in  this  section. 

THiL  Tjo  CA3:::3  of  TH2  DiFFi::i:]'TiAL  E;:uATion 

Equation  {?)  it^  the  general  e-^uation  v/hich  repre- 
i>ent3  the  systei.i  of  forces  wliich  o  xe  s-.iii..aed  at  the  valve 
stem.   Two  solutions  to  the  equation  were  obtained,  Case 
I,  in  which  the  K,  P,  tera  and  the  tank"  pressure  during  the 
run  v/ere  kept  constant,  and  Case  II,  in  v/5-iich  the  K.P, 
term  varied  with  the  valve  position  as  sho\m  in  Plate  II, 
and  the  tank  pressure  decreased  \d.th  tiiie  during;  the  run, 
Vfliile  Case  I  does  not  accurately  describe  the  forces  act- 
ing on  the  control  valve,  it  is  simple,  and  its  solution 
gave  the  authors  experience  in  operating  the  relatively 
complicated  computer. 

The  actual  com outer  circuit  for  Case  I  is  shown 
on  Plate  IV.   It  v;ill  be  noted  thr.t  the  Invariant  K,  P, 
terra  is  simulated  by  a  constant  voltage  fed  to  the  sura- 
ming  point  through  a  one  megohm  resistor.   The  initial 
conditions  for  the  equation  v/ere  set  into  the  apparatus 
as  follov«rs:   prior  to  zero  time,  sudtches  A,  B  and  C 
v/ere  closed,  so  that  the  condensers  across  amplifiers 
3  and  S  were  shorted,  and  could  have  no  initial  charge. 


It  is  evident    that  this   condition  must  prevail   since 
both  velocity   and    integral    control  iiiust  be   zero   at 
zero  ti;ae.      TbJ.'-   condenser  across  aiQplifier  4  v/as 
initially   cxiarged   up  to   the   volta^^^j   of   battery   D, 
which   hu3  the  effect  of  applying:  aii  initial  value 
of  Pg   at   tiiue    equal    to   zero.     Huns  were   started   by 

closing:   the    three  switches   siiaultaneously . 

It   was   aasuiiied  that   in  the  operation  of    the 
wind  tunnel    the    control   system  will   be   in  operation 
before   the    quic,k--.n;)enin^;  valve   is  opened.      Thus  when 
a  desired  pressure   input    is   introduced  into  the   con- 
trolle'r  vdth   the   nuick- opening  valve   closed,    the  con- 
troller will  oiX;n  the   valve   to  full  open  position   in 
an    effort   to   increase    the   settling  chaiaber   pressure, 
which   of   coui'se  vdll   remain  at    zero   pounds   per   squcu-e 
inch  ga^e.      It  will   be  necessary    to    tceep   the    integral 
control  from  operating/   prior   to  tJie  openin^j  of  the 
quick- opening  valve.     Failure    to   do    this  would  build 
up  a  lar^e   integral    control   force.      When   the  quick- 
opening  valve    is   opened,    the  settling  chaiaber   pres- 
sure reaches    its  jiaximui'ii  value   in  a  short    tine.      In 
this   problGxL,    it   is   aEsuraec.   ti-)at    tiie    tlru-  required   for 
li.o   openiijg  of   tne  quick-openiiie.   v^^lve   ana  lor   tiie 
settlint^  chamber   pressure    to    rise   to    its  mcixiauia  value 
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is   zero.      This   iissuinoticn   iri,)t-ses    a  mora    severe   oon- 
dition  on    ti.t;   u  ..'iiuioii^jr    liSi    .aotaolx^     tv.,i:jt.s,      Tl'ie 
Cictual  probleiP.,  Jiov/ever,    is   sir-iplified  tlu'ou^ii  uso 
of   til  is  ass arapt  ion . 

In  Case   I,    at.   ^i-^ix   in  '•' :      rt^s   7   turouch  12, 

tiie    initial    value   of  the  voltfc^e   of  the  term/p^xx  500 

InJ/ 

Is  plus   30 1   v^iiich  correiy^onds   to   an   initial   vtulue  of 

p  of  35.8  psi^  when  n„  is  594  psig/ft  for  a  supply 

tsnk  Dressure  (p.)  of  93  osir;.  TUis  vulue  for  in- 

t  ' 

Itl.ai  pressui-e   upset    is   aocut    the    ...laxijaach  tiiat  would 

be    encountered   under   actual   operating  conditions  of 

the  v./ind   tunnel.      This    is   accounted  for   by   the   fact 

thc;.t   the    settlins;^  chainber   pressure  iaust   be   above   a 

.Tiinimuiti  of    52  psij;^  for    a   supersonic  shock-free   flov/ 

through   the  nozzle    test   section  when  the    "dlffuser" 

is  in  place. 

On  ridte   7  is  the   circuit  dia£:;raia  for    the 

coi:ipat€;r   for   Case   II.      Equation   (7)    is  inodified   for 

this    case   to    take   into    accvUi'f-.    ^'  '•   ^'■ariation  of  K  ? 

4   d 

vdth   the   control   valve  position.     Fro.a  Plate   11^ 

K,P,  =  11.86-  max  '^' 
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From  e qua t lens    (2)    and    (3), 

CoJibiniii-;   equations   (8)    and    (9): 

K,Pd  =  11.88  -  H2(-^J-    IM2-E    +  I4Z-^^  (.10) 

Tliiis  equation    (y)    is  modified  for  Case   II  as   f-llov/c 

^ji^yi  .11.88  -  ,Hz(|-^)  -  i.iE  -  'H.  f;  =  0 


To  sirp-uiutt!  reduction  of   tiie    tank  pressure 

during   a  run,    all    the   terms   ttiat   contain  the    coefficient 

n2  must  be  varied  with    time.     On  Plate  VI  are   plots  of 

n   ,   n     ana  ^, /n     versus    supply   tank  pressure  which  v/ere 

derived   froiQ  Plate   I,      On  Plate   VII- are  found  plots  of 

p  /  p       for    the    suTinly   tank  versus  time  and      (T*     versus 
t       t^ 

time   for    two   different    <r'^  ,   with    cr'  defined  as   the 
ratio  of   initial   stora{;e   tank  pressure   to    settling   cham- 
ber  pressure,    or(p^    /p)    ,       .     It   can   be   s^'--'-.    ^\'-^^-  -^r-^'^ 

t\        abs 
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equals  1.25,    the   len,':;th   of  the   run  is  3.5   seconds. 
This    appears   to   be    about    as   sevi^re  a    ctiae   of   pres- 
sure drop   as  would  be  encountered   in  actual    operat- 
ion of   the  wind   tunnel,      ^^ssumin^  an  initial    tank 
pressure  of   93  psic;  «nd    a     o-'    of  1.25,    it  follows 
that    the  desired   settling  chamber   pressure   is  71.5 
psi{;,    and   the    tank  pressure  at   the    end   of    the   run 
would   be  about   75   psi^. 

^  voltage   to   give    an  initial   p^   of  21,5 

e 

psi{/,,   which   corresponds    to   a  desired   settling   cham- 
ber   prestiui^e  of   71.5  psig  was  introduced  into   the 
co/nputer  circuit  by   means  of.  battery  D  as   indicated 
on  Plate  V,      b^vitches   A,   B   aiid  C   operate    in   the    same 
iMnner  as   described   for  Case   I. 

Referring  then    to  Plate  VI,    if    tne   tank  pres- 
sure will  drop  from  93   to   75    osij  in  3o5  seconds,    then 
n2  must  decrease  froia  594  to  512  psi£;/foot   and  n^ /^p 
must  decrease   from  0.05^7  to  0.049  feet  durint'    the 
same  time   interval.      This  vai'iation  was  assoiaed   to  be 
linear  with   tiii^e.     From  Plate   VII   it  can  be   seen   that 

for    a   small    <r*  ,    the    variation  of  p,/p.      v/ith  time 

t     t  < 

can   be   approximated   hy  a   straiglit  line,    so  that   the 
above  as  sunpti  on  is  ne  arly    cor  rec  t , 

Referring;   to   equations (6)    and    (11)   and   to 
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Plate   V,    it  is   noted    thr^it  tUe    terrn  n,,    is  cr.nf. o  l.iu--.-^    .-js 

/^ 

u   direct  lU'-OTjor tionalitv   rixtvjr    iu   the.   teras  i:   ,    1: 

1*      2 

and  k-.    of    ttie  differential    equation,      Tiius,    by  vary- 
iric:   tile   voltat'e   output  of   eu'iplifier    6   in  direct   pro- 
portion to   ti>-;    VcUue    of   iip,    trie  variations  of   ICt  ,    k^ 
and  k-    that  were  caused   by    tiie    dropping  of  taiik  pros- 
sure  were    simulated. 

jn-lso  in  e^iuation   (11),    the  tc-ivij:i   142p/n2 
and   142n-/np  were  vaiied  v/ith    tiiae.      In  this  case, 
p  -.V6S  taken  as   71.5  psig  as  was  previously  derived. 
Then   at   zero   tLae,    triu    two    t<^rms  were  -17.1   i-Ui'i  3.3 
volts  repectively,    and   at  zero  plus  3.5   seconds, 
they  were  -19.9   and  7.0  volts. 

PotentiOiuetcrs  wer<=   eLiployea  as   in  Plate  V 
to  vary  simultaneously   the   volta/:es  described    in  the 
ubove  two   paragraphs  linearly  with    tiirie  .      The    three 
potentiome  ttrr  s  were    arranG*--^   iii   tcinde.a  on  a    siii^^le 
shaft,   or   ''ganged",    so   that   such   control   over  the 
effect   of  tank  pressui'e   variation  could   be    simulated » 

In   addition  to    ctppearin(-   as  a    coefficient 
in    several    of   the    terms   of   the   differential   equation, 
the  factor   n2   also  appeal's   as   a  part  of  the   variable 
^e     2   ^^^   ^^^  derivative  and    integral    terms.     No 
physical  operation  during  the    runs  was  employed   to 
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vary  the   factor  112   in   this  variable.      Its  linear  var- 
iation v/ith  tilde  v;as   taken    into   account  in  evaluating; 
the  or di nates  of   the   vaj^ious  functions   of    tht;   variable 
froiii  the  recording   tape. 

•    RLSU1,T3   ATD  DISCUSSION 

The    tiue  variation  of   the   various  pertinent 
terras  of   the   differential    equation,    as  obtained   by   use 
of   the  computer,  is  snown   in  Figures  7   tnrou^h  12  for 
Case   I,    and   in  Fit:ures  13    thi"OU^h  17   for  Case   IIo 

Values  of  k   ,    k     and    k^  were   chosen  experi- 
.Tientally   sr.ch   th  it    tne   pressui'e   error  decreased   to   a 
ne^licible    amount   in  the    shortc3St  possible  titne  with- 
out   exceeding    the    allowable    forces   tiiat    the   motor    i." 
capable   of   trans.Tiittini^   to  ti.c    valve   steia.      The  values 
of    the   resistances  l/k,  ,   l/iC2  ^nd  l/k^   that   were   used 
■*•'•    ■*•'':•    solution  of  the    dif  fe;;;ential    equation  for  both 
cases   ai'e  found    in  Table    II.    /J. so  in   this   table   are 
the    corresponding?;  values  of   the  coef ''icients  k    ,    k^^ 
ana  Ko  .      These   coefficients,    when  xaultiplied  by   the 

instantaneous  values  of  d    ,    p.    <'ind    |  p  dt   respectivelv  , 

^e        e  Jo    e 

describe    the    forces   amlied   to    the    valve    stem  at   any 
instaiit  of   time   by    tne    acrivaioive.  control ,    proportional 
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O'"^nt"r»ol  and   integr'T    <"  vni.r:"'^    :^  '    t  w'^e'^-    ■"'f  ■^T^c^    •-<^-^•1|,■,-— 
iu  Chan  ism. 

No  attempt  v/as  made  to  design  and  build  a 
servoiaechanisii  f -^r^  tiiis   cont.rol   p^oble!   d')P  t-.,-,   t':^- 
time   limitationc.      Hov/ever,    steps  were  taken   to    Su-' 
'.ain  liiaterials  and   equipment  tov;ard  this   end.      A 
control    v--lve   si-'i-^  "^/^r  t;'>  tVi,-it,    nrr-.^pntly  installed, 
':s   oeen  in  Figure   i,   and  a   ocrvoaotor,   the   character- 
istic 3  of  vfhich  are   shown  on  Plate   IX,   have  been 
procured.      It   is   planned   th  t   the   project  i/ill  be 
carried  to   completion  usinK;  tr^ese   iteius  and  the 
results  of  this   study. 

^'•e.viltn   of  the   G-jsn    I    s^olution  as   obtained 
from  the  various   .;rush  recorder  x,apec>   are   as  f  ollov/s : 

"''ifure  7(a)  :  1:; xi:.ia:i  requir^sd   control   force:  13    lbs. 

required    .steady  state  control   force:    5.3    lbs. 

(b,  p^_  initial:  35.7  psi. 

p^  at  t  =  l   second:  3.0  psi. 

6  ... 

p  j.  at  t  =  1.5   seconds:  0.0  psi, 

p     at  t=2    seconds:  1.2   ^-^si. 

^e 

p     at   t  =  2.2   seconds  0.0  psi. 

"igure   3(a)  Coulomb  friction  varies  from  +  1.23   lbs. 

to -1,23    lbs.   dr:nGnd5.nr  on  the  direction 
of  motion   of  i:a&  valve.        '"owever,   at 
about  t  =   0.7  second  j    amplifier  7 


appai-eiitly  failed   to    saturate,    so   this 

term  is  in  error   ptst   that   time.      It 

v/ill   be   sJiovvn  later   that    the   coulOiub 

friction  effect  is  relatively   aai-;!!  ^/h^^n 

coiTipared   to   the  othwr    effects, 

(b)  This  is   the    3b.ne    trace   as  Figure   7(b)  , 

and.   was    Ltjuxua  ^    i  ^...     '^viLie    cOiUparisons, 

Fi^rure   9{a)  MaxLiiuiii  re-xuired   control  force:  8     Iba. 

Required  steady    state   control 

force:  o.plbs. 

This  plot  and  Ficure   9(b)   v;ere   included 

only   as  a  inetter   of    interest.     They  re- 

Dresent    the   cor.dition   in  which  the   K  P 

4  d 

tennis   positive.      This  v.ould  corres- 
pond  to  a   lift   force   on   the   control   vulve 
ixi   tiie   positive   or    dovm  directioi^  ,k-i.<su. 
vras  opposite    to    the  lift   force   as  measur- 
ed, 
(b)  p     initial:  3^.7  pai. 

6 

p     at   t=l  second:  2,4    osi. 

p     a^   t2l.5   seconds:  C,9   t'sI. 

e 

Fi-ure   10  /.  > 

{^)  (p   /:  .,)  iiiaz  =  (i)  iiax:  1.1   ft/sec  =13.2    inch/sec 

[p  )  max  at  g=  0,12  second:         6^3  psi/sec, 

p  at  t=l  second:  negligible 
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(b)  This  is   the   saoie   trace  as  Figure   7(t))  , 

and   was   included  for  time  coiiiptirison. 

2 
Figure  11(a)  I'^q^^z)    ^J^^  =  (^)    ^^^-  22  ft/sec    . 

Refer    to   equtition    (5)    for    the   above  re- 
lationship.     The   iaaximuia  acceleration 
occurs   at   zero  time.      It  will  be  noted 
that    the   effect   of    ttie    couloiub   friction 
decreases    the   acceleration  slightly  when 
the    velocity    chan^^es   direction, 
(b)  This  is    also    the    same   trace  as  Figure   7(b), 

and  was   included   for   ti^ie    coinx^arison. 
FiGure  12(a)  This   is    the   saiae    trace   as  Fis^^re   7(a),    and 

v;as    included   for   tiine    coraparisono 


Jo 


(b)  Maxiaiuxi         p  dt  at   t  =   1,5   seconds:    8.3  lbs. 


[p  dt   steady    state:  6^8  IbSo 

It   is  noted  that   in   the   steady    state   the 

only  control   force   tictin^  is  that  v^'hich  is 

'due  to   integral  control.     Comparing  the 

steady   state   integral    control   force   of   6. Bibs. 

with   tiie    steady   state   control   force  of   5*3 

lbs.   from  Figure  7(a),    ^  discrepancy  of 

1.5  lbs.   is   found.      This  disagreement    is 

Dossibly  due    to   Uie    f^ct   that   p     and   p 

■^e  e 

differ   from  zero  by  a  small  amount   at   the 
time    selected  for   steady    state  comparison. 
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The  Brusii  recorder  tap^s  on  wiiich  the  solut- 
ion to  Case  II  v/as  recorded  -ave  inrorjiL-tion  as  belu\\ . 
It  should  be  ruJie.abored  that  ordinatea  that  contain 
the  factor  n^   in  the  dependent  variably  idust  be  correct- 
ed for  the  variation  of  n^  with  time.  Plate  7III,  which 
is  derived  from  Plate  VI,  is  an  aid  in  determining  the 
value  of  D.2   at  any  time  t. 

Figure  13(a)    Maximum  required  control  force:     i^u   lbs. 

Required  steady  state  control  force :13  lbs. 
(b)    p^  initial:  21.5  psi. 

p^max  at  t=  0.75  second:  64,5  psi. 

p^  at  t  =  l  second:  3,4  psi. 

Pg  at  t  =  1.2  seconds:  ne^-li^-ible 

Figure  14(a)    The  coulomb  friction  force  varies  in  a  manner 

similar  to  th^t  of  Figure  8(a).   ijaplifier 
7  apparently  failed  to  saturate  at  about 
t=0,6  second  in  this  case, 
(b)   .  This  trace  is  identical  to  Figure  13(b), 
and  was  included  for  time  comparison. 
Ficure  15(a)    (p^/n2)max  =  (i)max:    I.3  ft/sec  =15.6  in. /sec. 

(p  )max  at  t  =  0.09  second;.       770  psi/sec. 
Pg  at  t  =  0.6  second:  negligible, 

(b)    This  trace  is  identical  to  Figure  13(b),  and 
v/as  included  for  time  comparison. 
Ficare  16(a)    (p  /n_)maz  =  (x)nax  at  t  =  0.0  second:  43  ft/sec'^ 


Tho   small   effect   of  coulomb  friction   is 
noted   i.'.'.   u.is   trace   ulso. 
(b)  This   trace   is   identical   to  Figure  13(b), 

end  was   Included   for   time   com.oarison. 
Figure   17(b)  MaxiinOxfr    [p  dt    lo  c/c    t=   1.."    ^^<:-o^.'^d:    8.15   lbs. 

ELEGTRCKIG   Di:'yLRLNTI.iL   I^hd^YZlLR  DISCUSSION 

Co.nplicatfcd   circuits   should  not   be  attezaptbc  by 
those  unfamiliar   v-'it  i   coLiouter   oiperution.      It   is  recom- 
'mended   that   -n  simplified  vc-rsion  of   the   probleiti  be   solved 
first,    so   that   a   der^ree     ""   f suniliarit^   with   the   ooiaplej.- 
ities   of   the   computer  iniciht  be   obtained   by   the   oper&tor 
before   proceeding   to  iaore  difficult   circuits. 

The   usual  matheittatical   laethods   for   solving 
second  order  differential   equations  such  as  are   contciined 
In   tbJ3  report  v;ould  prove   to  be   extremely   laborious  and 
tinie   consuffiint:.      Cn  the   other   hand,    the   computer   seems 
to  be   an  ideal  means   for  obtaining  such  solutions.      The 
-Tiost  difficult   part   of  ,the  procedure   is   the    settint^  up 
'?^   ^^e  equipiaent   to  i3uplicate   the  equation.      Tlien   the 
values   of  the   constants   are    obtained  by  a   trial   and 
error  process,   notinij;   the  response  of  the  oscillograph. 
v;hen   3   3atisf{^otory   ^^'^^  v-.n.;^::    h;is  >.H^--n   obtain^^f] ,    the 
constants  which   caused   this   response   are   knovvn,    and   the 
problem  is  solved.      Tiie   constants   thus  determined  can  be 
usee*   to  design  the  actual   control    system. 


Resistors  and  shorting  bars  were  mounted  on 
General  Radio  type  274  -  M  double  plugs  whicJa  fitted  into 
ttie  terminals  of  the  amplifiers.  Erratic  results  from 
time  to  time  were  finally  traced  to  the  looseness  of  the 
mechanical  connections  between  tiie  resistors  and  plu^s. 
Time  can  be  saved  by  checkin,-  these  connections  before 
insertlni^  the  plut'S  into  the  amplifier  terminals, 

Tne  amplifiers  used  ore  of  an  inexpensive  type 
constructed  several  years  ago  and  must  be  individually 
balanced  for  zero  output,   wsdlcv/iiii;^  the  amplifier  tubes 
to  heat  up  for  a  long  period  of  time  has  a  favorable 
effect  on  maintaining  balance.   It  i^  recommended  that 
the  tubes  be  allowed  to  operate  overnight  if  possible, 
Rco^n.   temperature  and  hoiaidity  also  effect  balance,  and 
where  tneir  variation  is  not  closely  controlled,  frequent 
balance  checking  is  necessary,   Cf  course,  a  well  regulat- 
ea  power  supply  for  the  amplifiers  is  also  necessary  for 
F.ood  balance. 


CONCLUSIONS 

The  objective  of  this  study  was  to  determine 
the  characteristics  of  such  a  control  systtjm  cis  vvoulu 
reduce  the  pressure ''error  in  the  settling  chamber  of 
a  particular  blowdown  v/ind  tunnel  to  within  plus  or 
minus  one  pound  pwr  square  inch  within  one  second.   By 
referrint^  to  Figure  13  it  may  be  seen  that  the  solut- 
ion obtained  in  this  study  can  be  considered  as  in 
satisfijctory  agreement  v.'ith  the  above  stated  object- 
ive. Tne  values  of  the  constants  associated  v/ith  the 
solution  are  found  in  Table  II, 

It  is  recoruinended  that  the  design  and  fab- 
rication of  a  servome Chan ism  that  will  exhibit  the 
characteristics  as  set  forth  in  Table  II  be  accomp- 
lished, so  that  satisfactory  pressure  control  in  tii.e 
subject  wind  tunnel  nxa;/   be  obtained. 
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coMPumnoNS 


DRIVE  COMPUTATIONS 


Data: 

Weight  of  V'-dve    and  steipj  10o38   lbs. 

Coulomb  friction  force i  1.23   lbs. 

"Sticktion"  force:  approx.1.00   lbs. 

Valve   aerodynamic  lift  See  Plate   II. 

Assumed  config-arations  See   Fig.   2. 

Acceleration  forces      Assuming  a  maximum  acceleration  of  6   ir./aeo    « 
Force  sr    m  x  a  ss  ^^^°|?^(0o5)     »    0.254    lbs. 

Ibtal  Force  1     Neglecting  Lift  Force   other   than  the   6   lbs.    for  v.hioh 

the  valve    is   compensated, 

Fr.=  F    +?.-fF;,    S3    0.25  -i- 1.23  i- 1.00   »\4B    lbs. 
^        A  V  b 

Work/sec  required  of  motors      2.48(0.75)  -^    0.31   ft.    lbs. 

(I2)i.t)  seo. 

Required  Horsepower:  0.000564  HFo 

Motor   Data:      Diehl,   Navy   type   CDA   ?,11.')52,   Max.    output  =  5,2  watts   or 

0,00693  HP  ®  5.5  oz.    in.    torque    and   2000   rpmo 

If  motor    is   to  aocelerate   between  zero  and  2000  rpm  in  ^  sec^rid, 

the    torque  'would  vary   from  5.46    to   3.5   oz.    inches.       This  would  be 

equivalent   t'^   operating  at  4.7   oz.    inches   and   lOOO  rpm  constantly 

for   one   half   gee end. 

Work/sec  =(4.7) (1000) (2  x  3.14)     -     2.56   ft.    lbs. 
(1-.  }(16)    (60)  ~"'6ea 

This  compares   to  0.31   ft.    lbs, /see  required  by   the   2.48   lb.    load. 
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Ihus   the  motor  would  deliver   a  force   of  20.5   lbs.    at  the  valve 
stem.  X       «    2,48  x    :»    20.5    lbs. 

TTo  6  "  o73T    ^ 
In  order   to  exert  a  20.5   lb.    force   on   the  valve    stem,    the 

average   torque   of  4.7   oz.    inches  must  act  through   a  gear    train: 

(4.7)  X     s.    2:^.5       ,        X  »  6y.7      gear   ratio. 
TT6)  ^ 

After  a  solution  to  Casel  was  found,  greater  control  forces 

were  required  to  give  a  larger  margin  of  available  control  motor 

force.   Ihe  Diehl,  Navy  typ«  CDA  211426  servomotor  was  available, 

which  would  operate  between  8.4  and  11,5  oz.  inches  over  the  same 
rpm  range.  'Ihe   equivalent  constant  torque  in  this  case  is  1"'.4 

oz.  inches,  and  the  motor  would  deliver  an  average  force  of 

56.7  pounds  to  the  valve  stem.   Ihe  gear  ratio  necessary  for 

this  motor  would  be  87  to  1.     (b6.7)(16)  ^  87,0  gear  ratio. 

TTcuTJ 
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I-rinoiijl'i  ijir-cnaions   of 
University   of  Iiichi^^an  Blov.'  Dovm 
V/ind  Tunnel 


Volaiue   of   storii^fej    tanK : 

100         ft.'' 

KaiLi;na:a  supply  line  prossurt  : 

95          psiB. 

^'^::■l^    throfct   S'3Gti".     •■---r,  • 

0.28  ino^ 

x\o-i.7.1e   test   section  area: 

1.00  in.'^ 

A"!l    :^la!iibinc^  betweer.    i' -■  ■ 

...nd  no'zy.le   is 

3   iiicii  pipe   and  fixtuios. 

'iiBL.j  I 


-i"d~  ' 


Control  Force  Constants 


ohms 

OllillS 

oiiiiis 

k 
1 

k 
2 

ft 

S 

^1 

i   sec- 
psi 

K  , 

psi 

K3 

ft     360 

pDi  sec 

Case   I       0.3 

0.008 

0.003 

3.33 

125 

333 

0.0056 

0.21 

0.56 

C:a£;    II     0.15 

0.003 

0.001 

6.67 

333 

1000 

0.0112 

0.56 

1.68 

TxiBIi:;  II 


-J/- 


Settling  Chiiiiber  rresjiu-c   vs.    Valve  Position 


X 

inches 

0.125 

0.375 

0.625 

0.375 

1.125 

lo375 

Pm 

93  psig 

Al.O 

53.5 

66.0 

76o0 

79o0 

80.0 

\ 

BS  psi^ 

36.8 

47o8 

59o2 

6?o4 

72.3 

73  o2 

\ 

76  psi- 

33.3 

43.0 

54.5 

60.0 

63  o3 

63  o5 

Pt. 

1 

60  psig 

24.3 

32.5 

42.5 

43. a 

49«5 

51.0 

P  = 

psis 

TiiBIiC   III 
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Lift  Force  vs  Valve;  rosition 


X 

in. 

Lift 
Ibo. 

0.5625 

10.19 

0.6875 

12.69 

0.750 

13.50 

0,675 

14.50 

loOOO 

17. CO 

1.125 

17.175 

-   T  7  C 

20,00 

TABLE  IV 


-33  - 


Detemination  of  Coulomb 
Friction  Force  and  '..'eight 


1" 

Up  I  lot  ion 

)oyrn  LTotion 

lbs. 

lbs. 

11.500 

9.000 

11.6S7 

9.31? 

ii.au 

9.063 

11.406 

9.313 

11.656 

9.C63 

11.625 

9.ia^ 

11.563 

9.156 

'vj:      11.612=V;4 

? 

Avg : 

9.156»W-.F 

v; 

as 

10.3^  lbr3. 

F 

= 

1.23 

lb-:. 

rp    •<   ~T   T?       TT 


—5-^- 


/=7^  J. 


-  3-^- 


471  . 

JQ     Reduction  6«*r,  697:/ 


/ 


/ 


Sketch    of    ProjDoseol     Dirit^e    MechAnIsm 

F.G-.2 


Sketch  Xncltcatin9  Method  Used  to  Measura 
Coulomb  Fyiction  ^^nd  Vki<^hi  of  Vzlve 

Ro.  3 


v°^  "^'T  *>^ 


rr^ 


..  <N:n^  . 


Sketch  Xndicdtino  Method  Used  to  Measure 
"Lift"  Force  of  l/alve 


valve  Closed 

1. 

5" 

+ 

All  displacements  and  forces 

are  measured  bositivc  ^rom 

the  closed  to  the  of>ened  |3osition 

[ 

\/a  1  &/./%         ^^tf%^%  ^^^k.^ 

vaive   upeneo 

Coordinate    Si^stem 

R&.5 


—^7-^ 


^'J-  ^ 


—  J?-^— 


Fig.   7(a)  I     Control   Force.      lAi*  0»3J    ^Ag*  0»008j   1A5*  0«008, 

all   In  megohms. 
Fig.  7(b)  I     £e  x  500»     K.Pda-6  volts.     Constants  same  as  7(a). 


Fig.   8(a):      Coulomb  Friction  Force.   Constants   seme   as   Fig.  7(a). 
Fig.   8(b)  I     Same   as   Fig.   7(b). 


—  ^^- 


Fig.  9(a) t     Control  Foroe.  Constants  same  as  Fig.  7(a). 

Fig.   9(b)  t     2fi  X  500.     K.Pd«-*>6  Tolts.     Constants  same  as  7i^a). 


Fig.   10(a)  I     2b*  Constants   sane  as  in  Fig.  7(a). 
Fig.   10(b)  t     Saiae  as  Fig.  7(b). 


-7^<0    — 


Fig.   11(a)  J     Be,     Constants   same   as  in  Fig.   7(a). 
Fig.   11(b):      Same  as  Fig,   7(b). 


Fig.    12(a)  J      Same   as    FIr.    7(a). 

Fig.    12(b) »      Integral   Force  x  500.   Constants   same   as  7(a). 


—    ^/~ 


Fig.   13(a) »     Control  Force,     l/kj^^     O.I61   lAg     O.OOSi   l/k^     0,001, 

all   in  megohms* 
Fig.   13(b) I     §®  X  500.     Constants  same  as  Fig.  IS (a). 


Fig.   14(a)x     Coulomb  Friction  Force.   Constants   same  as  Fig.   13(a). 
Fig.    14(b)»      Seme   as   Fig.    13(b). 


-  ^2 


Fig.  15(a)  I  £».  Constants  sane  as  13(a).   (t)i  Same  as  13(b). 
^2 


Fig.  16(a)«  20.     Constants  same  as  13(a>.   (b)i  Saine  as  13(b). 

»2,.  .......  — -XT-T 


Fig.  17:  \-^   d*  ^   500.   Constants  sane  as  Fig.  13(a), 
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